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Abstract: Supersymmetric QED hydrogen- like bound states are remarkably similar 
to non-supersymmetric hydrogen, including an accidental degeneracy of the fine struc- 
ture and which is broken by the Lamb shift. This article classifies the states, calculates 
the leading order spectrum, and illustrates the results in several limits. The relation 
to other non-relativist ic bound states is explored. 
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1. Introduction 

Supersymmetric bound states provide a laboratory for studying dynamics in supersym- 
metric theories. Bound states hke hydrogen provide a framework for understanding the 
quahtative dynamics of QCD mesons, a supersymmetric version of QED can provide a 
quahtative picture for the symmetries and states of superQCD mesons. Furthermore, 
recent interest in dark matter as a composite state, leads to asking how supersymmetry 
acts upon these composite states [4-7]. 

This article calculates the leading order corrections to a hydrogen-like atoms in 
an exactly supersymmetric version of QED. Much of the degeneracy is broken by the 
fine structure and a seminal calculation was performed in [1] for positronium, see [2] 
for an A/" = 2 version of positronium. Supersymmetric hydrogen is a similar except 
for the absence of annihilation diagrams, see [3] for an independent calculation. In 
the heavy proton mass limit, the supersymmetric interactions of the theory become 
irrelevant operators, suppressed by powers of the proton mass like the magnetic moment 
operator in QED and the fine structure is identical to the non-supersymmetric theory. 
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This article finds that fine structure spectrum of supersymmetric spectrum of hydrogen 
has an accidental degeneracy which is exactly analogous to the accidental degeneracy 
of the I = and I = 1 levels of the n = 2,j = | state of hydrogen. The supersymmetric 
version of the Lamb shift lifts the residual degeneracy and this article computes the 
logarithmically enhanced breaking. 

The organization of the paper is as follows. In Sec. 2 the general framework for 
supersymmetric QED is set up. Particular emphasis is placed upon the symmetries 
and states of the theory. The non-relativistic limit is taken in Sec. 2.2 and in Sec. 
2.3 the equivalent of the heavy proton limit is taken. The accidental degeneracy is 
demonstrated in this limit and the organization of the states is illustrated. Sec. 3 
calculates the fine structure of supersymmetric QED and the results are similar to [3]. 
The Lamb shift is computed in Sec. 3.1 and it is shown to break the residual degeneracy. 
Non-Abelian Coulombic bound states are discussed in several different circumstances 
in Sec. 4 and the properties that carry over from hydrogen are discussed. Finally 
applications are discussed in Sec. 5. 



2. Supersymmetric QED 

The field content for supersymmetric QED consists of four chiral superfield: E, E'^, 
P, P'^ with charges: —1, +1, +1, —1 under an Abelian vector superfield, V. The 
Lagrangian for supersymmetric QED is given by 

/:SusyQED = jd''eK + jdPeW+ h.C. (2.1) 

where K is the Kahler potential and is given by 

K = ^^eMq^gV)^ (2.2) 

where is the charge of the chiral superfield, $, and g is the gauge couphng of the 
U{1) vector superfield. The superpotential is given by 

W = rrieEE" + rUpPP" + ^W^W^ (2.3) 

where rrie and rrip are electron and proton masses, respectively and Wq, is the super- 
symmetric gauge field strength. The Lagrangian Eq. 2.1 is expanded and is given 
by 

-^Susy QED = C,Km.+ Mass + -^Susy Int. (2.4) 
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where £Kin.+ Mass IS the Lagrangian for the kinetic and mass terms of the theory and 
-^Susy Int. contain the supersymmetric interactions. These are 

-^^^11.+ Mass = eipe + + ipe" + meiee" + h.c. ) 

+ |D^e|2 + |Z}^gf + m^(|g|2 + |gf )■ 

+ terms for P and 

-^-F'^^ + M^X (2.5) 

where e and e denote the Weyl fermion and complex scalar from the superfield E and 
Ffj^i, is the gauge field strength of the vector field and A is its respective gaugino. The 
gauge interactions are contained inside the covariant derivatives, D^. The supersym- 
metric interactions are 



>Csusy Int. = iV2g\ [ee^ H ) + h.c. 



-^(|eT-|ef -Ij^r + rn'. (2.6) 
All of the interactions of the theory are proportional to the gauge fine structure constant 

and throughout this article, weak coupling is assumed, a <^ 1; though not necessarily 
a = ctEM = 1/137. 

This article considers the bound states of an electron (or its respective superpart- 
ners) and a proton (and its respective superpartners). Throughout this article, the 
relative masses will satisfy rrie < rup. The remaining portion of this section contains 
a description of the symmetries and the states of the theory in Sec. 2.1, the non- 
relativistic limit in Sec. 2.2, and the heavy proton limit in Sec. 2.3. 

2.1 Symmetries and States 

The Lagrangian in Eq. 2.1 has all of the symmetries of QED: rotational invariance and 
a Z2 spatial parity. In addition to parity, the Lagrangian is supplemented with a U{1)r 
symmetry. Parity plays a major role in the organization of the spectrum and it acts as 

(-l)^a;'^ r 4-^ V ^-V E^E' P ^ P^ (2.8) 

This implies the supersymmetric derivatives, and field strength, Wq, = 'D^T'aV 
behaves as 

^ (2.9) 
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Most significantly, the U{1)r symmetry does not commute with Z2 parity and are com- 
bined into an 0(2)ij symmetry. The 0(2)ij symmetry is important because it has two 
dimensional irreducible representations. The fine structure spectrum of supersymmet- 
ric hydrogen will contain many two-fold degeneracies; however, half arise because the 
states are doublets of the 0{2)ji symmetry. 

In addition to the space-time symmetries, the constituent particles are also charged 
under a global f/(l)e x U{l)p flavor symmetry and the states of hydrogen are charged 
under the [/(l)e+p diagonal subgroup. The spectrum of supersymmetric hydrogen is 
doubled relative to supersymmetric positronium because hydrogen is charged under 

The leading order interaction between the supersymmetric electrons and protons 
arise from the Coulomb interaction and the states organize into quantum mechanical 
hydrogen spectrum, with state \nlm), given by 

i^nimix) = {x\nlm) = Rni{r)YimiO,(f)) n>0 l<n \m\ < I (2.10) 

where Rni{r) are the associated Laguerre polynomials and Yim{0,(j)) are spherical har- 
monics with angular momentum /. In addition to the bound state spectrum, there are 
also the continuum states that play an important role in guaranteeing a supersymmet- 
ric answer for the fine structure spectrum [1]. The energy spectrum at leading order is 
given by 

£r- = (2.11) 



where /i is the reduced mass 



fi-^=m~^ + m~\ (2.12) 



The fold degeneracy at each level is a result of an accidental enhanced 0(4) sym- 
metry of the non-relativistic Coulomb potential. This is broken relativistic and super- 
symmetric effects at 0{a'^fi). In addition to the spatial wave functions, the spin and 
supersymmetric portions to the wave functions, \S) and is 

\S),^^^ = {E® i5^^)p=p.^+p,,/2 ® (P © P^^)p=p..-p,,/2. (2.13) 

Non-relativistic bound states are necessarily massive and therefore fill out massive N=l 
multiplets which are built from the Clifford vacua, Qj, of spin j. These have the field 
content of massless N=2 multiplets. Clifford vacua are useful for classifying the excited 
spectra in Sec. 3. This tensor product can be decomposed as 

\ y Massive charged vector multiplet 
\S) = { ^ ^ . (2.14) 

I J^2 Two massive charged hypermultiplets 
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where contains the states of a massive charged vector superfield and ^2 contains 
the states of two massive charged hypermultiplets, which is equivalent to jour chiral 
superfields. These states are charged under the C/(l)e+p symmetry of the theory unhke 



positronium. 


Exphcitly, the states of 'f and 


are 






massive complex vector 


a;_i_,a;_ 


r = < 




massive Dirac spinor 






X2,X2 


massive Dirac spinor 


C2, C2 






massive complex scalar 


^ ■CC7-|_, VJ 



massive complex scalars 
massive Dirac spinor 
massive Dirac spinor 
massive complex scalars 



(2.15) 



The 0(2) symmetry forces ro, x, and C into doublets. Parity acts as 



c. 



— V 'Vf, xi ^ X2 Ci 
on the fields with spins and transforms the scalars as 



(2.16) 



(2.17) 



The Coulomb potential is insensitive to the spin/supersymmetric state and therefore 
the states factorize at 0[a^^) as 



IV^) = \nlm) ® \S). 



(2.18) 



The tensor product between \lm) and \S) is performed and the states organize them- 
selves into massive supermultiplets. Using the notation VLj to denote a massive super- 
multiplet built from a Clifford vacuum with j, e.g. VLq is a hypermultiplet, the the 
tensor product between a state with orbital angular I and a multiplet built from is 
given by 



\lm) ®VLj = Vt 



The states of principle quantum number n are organized into states 



\i-3\^^\i-3\+i (2.19) 
nized into states 

(2.20) 

where there is a two-fold multiplicity of states for every VLj except j = n — \. The 
leading order superspin wave functions can be calculated by acting with the susy rais- 
ing operators on the Clifford vacua in the expression above. The states where the 
constituents arc m a, Y configuration arc mixed with each other because the Clifford 
vacuum has non-zero spin. There are two basis choices 



4.(10001) vs |0®(aL-Qi) 



(2.21) 
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where aj^ is the supersymmetric raising operator. The first expression in Eq. 2.21 
gives the proper organization of states. This leads to mixing between states where 
the constituents are in the and ^_ configurations. The relative admixtures can be 
computed with Clebsch-Gordon coefficients. 

The two-fold multiplicity of integer j arises from these states are doublets of 0(2)/j, 
whereas the degeneracy for half-integer j is accidental. For the ground state of super- 
symmetric hydrogen, there is no degeneracy of the spatial wave function and therefore 
the states are 



Pi = 2 X fio © fii (2.22) 
where ~ fii and ^ = 2 x • 
2.2 Non-relativistic Limit 

Hydrogen is a non-relativistic bound state and it is convenient to take the non-relativistic 
limit of the Lagrangian in Eq. 2.4 in order to gain intuition about the fine structure. 
The non-relativistic limit of fermions is found by organizing the Weyl fermions of Eq. 
2.4 into Dirac fermions 

*?'=(;.) *r=(;o). (2.23) 

After transforming from the Weyl basis where 70 is off diagonal to the Dirac basis where 
7o is diagonal, the non-relativistic limit is found by factoring out the fast exp{imt) of 
the wave function and integrating out the small components. For instance, the non- 
relativistic spinor for the electron, ipg, is related to the Dirac spinor by 



(2.24) 




and the leading order Schrodinger Lagrangian is 

- 

NR = Ipel'dti^e + Ipe- Ipe H • (2.25) 

zrue 

The non-relativistic limit for the scalars is found by first factoring out the exp(zmt) 
and then by rescaling the wave functions 

g = e*™^*— g'^ = e*™^*-^0e- (2.26) 
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The non-relativistic scalars have the same engineering dimension of | as non-relativistic 
fermions. The scalars also have a Schrodinger Lagrangian 

Ce NR = <Plldt(t>e + <t>\^^e + " " " • (2.27) 

The key aspect to non-relativistic limit of supersymmetric QED is that the super- 
symmetric interactions become higher dimension interactions, much like the magnetic 
dipole of QED. In the non-relativistic limit, the supersymmetric interactions of Eq. 2.6 
become 

Csnsy Int = fi-A ( ^^(?/^e0e + C0eO + " = (^p0p + V^p0pO ) + ^.C. 



J [^m' - \^eA') - ^m' - ) ■ (2.28) 

The leading supersymmetric corrections to the spectrum will always use one interac- 
tion involving an electron state and one involving a proton state and therefore all of 
these interactions are suppressed by powers of (memp)~5 for gaugino interactions and 
{mf,mp)~^ for D-term interactions. The next subsection considers the heavy proton 
limit as a simplified example to gain intuition about the full calculation. 

2.3 Heavy Constituent Limit 

This section considers the limit where the constituents become much heavier than the 
other 

Tfl 

— ^ ^ 11 = me. (2.29) 

nip 

In this limit the supersymmetric interactions do not contribute to the fine structure 
of the supersymmetric hydrogen atom. This is seen from Eq. 2.28 where all the 
interactions between an electron state and proton state are suppressed by factors of 
(memp)~", n > ^. Sec. 3 computes that fine structure and shows that n = 1 because 
gaugino exchange only contributes at second order in perturbation theory due to parity. 
There is no way of obtaining any positive powers of rrip in the numerator and this 
immediately implies that fine structure of supersymmetric hydrogen given by the purely 
QED corrections to the Coulomb potential. These are exactly solved for [8] and are 

E„r = = = me ^ + ^ r 2.30 

I / a ^2 2n2 2n' \An j + ' 
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where j = I for scalars and j = | / = / — ^, / + ^ for fermions. Sates with scalar 
electrons are not degenerate with the states with fermionic electrons. The difference 
arises from the Darwin term which doesn't exist for scalars as Zitterbewegung is a 
property of Dirac fermions. 

This spectrum immediately implies how the wave functions of the composite su- 
permultiplets are organized for mass eigenstates: the light state is either a scalar or a 
fermion and the different states of the supermultiplet are created by toggling the spin 
of the heavier state. For instance, the ground state for the supersymmetric hydrogen 
atom with an infinitely massive proton is given by 

=^2 = (^p © </>p © 0;c) ® (0e © CO y = ii^p ® (l^p ® <f)*pc) ^ (tpe) (2.31) 

and the mass splitting is given by 

E,n. - E,n, = Ef^ - eT = ^ (2-32) 



2 



where EnUj denotes the energy eigenvalues of the nth principle state built from the 
Clifford vacuum Q~ using the correct basis Eq. 2.21. Each of the Qi multiplets have 

■' 2 

the same energy eigenvalues, which is dictated by the 0(2)^ symmetry. 

The n = 2 level has an accident degeneracy which is exactly analogous to the 
/ = and / = 1 degeneracy of the Dirac at the n = 2 level. One of the Qi states 
arises from the electrons in an / = spatial configuration, while the other arises from 
the electrons in an / = 1 spatial configuration. These states have different parity and 
there is no symmetry that interchanges the two. In the normal hydrogen atom, this 
degeneracy is broken by the Lamb shift. In the supersymmetric hydrogen atom there 
is a possibility that the degeneracy is broken by supersymmetric interactions. In this 
case, for the supersymmetric interactions to break the degeneracy, there would be a 
supersymmetric hyperfine structure scaling as 



4 2 

am 



-^2^'=° ~ E^Qi^i ~ (2.33) 



^2n'f" -^2f7'= 

2 2 ''^p 



This article finds that this is not the case and that the degeneracy exists even with a 
finite proton mass and this accidental splitting is broken by the supersymmetric Lamb 
shift and the calculation is presented in Sec. 3.1. 

2.4 Finite Mass Wave Superspin Functions 

The superspin wave functions \S) are decomposed in terms of the elementary con- 
stituents to leading order by demanding that the multiplets f and J^2 close under su- 
persymmetry transformations. The supersymmetry transformations are approximated 
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on-shell, in the non-relativistic limit for the electron states as 



/^(Cae^a^e + CaCO, (2-34) 



where a^^ = Saa and = S'^^e^a- The transformation for protons are the similar. 

The superspin wave functions for a finite mass proton are computed in the prod- 
uct superfields: P E , P'^\ E'^\ P E'^^ P^^ E . After evaluating these product superfields 
on-shell, all of the independent degrees of freedom reside in the last superfield: P'^'^ E. 
The superspin configurations for and Y are calculated by acting with the super- 
space projection operators Vi,V2,'Pt that project on to chiral, anti-chiral, and vector 
superfields, respectively. Evaluating these superfields in the non-relativistic limit gives 



XuXt 




= Ce(j)pllJe - Se1pp(j)e 






= co(l)*pcipe - seipp(i)lc 




<?- = 






a;+ = 


(0p0:. + 0;.0e)/V2 






C2e{4>p<P*ec - (t)*pc(t)e)/ V2 + S20V'p^e 


Ci,cr 


V'ci = 


CeV'p^e + Se(t)p^e 




^C2 = 


cei)p(t)l. + se(l)*pcipe 




= 


(0p0e + 0;c0:c)/y2 




ZU^ = 





(2.35) 



(2.36) 



where the righthand equations label the on-shell degrees of freedom and the mixing 
angle is defined as 



tan 6 



rUe 

m. 



(2.37) 



p 



The ^ ^ limit of these superspin wave functions demonstrates the separation of the 
valence particles into the spin ^ electrons in Y and the spin electrons in Jif2 described 
in Sec. 2.3. 



3. Fine Structure in Super symmetric Hydrogen 

This section presents the fine structure of supersymmetric hydrogen. Therefore, the 
full analysis will not be presented in this article and [3] provides a detailed description. 
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Instead, a simplified calculation is presented that will illustrate the organization of the 
states and only uses first order perturbation theory. 

The standard method to calculate the energy spectrum to bound states starting 
from a field theoretic, effective Lagrangian is to use Breit's formulation of the effective 
Hamiltonian [9]. For instance, in QED the one photon exchange gives rise to a Born 
amplitude of 

= -v/(2Ei)(2E2)(2E3)(2E4) [wlwl UipuP2,q} ^lU'a] (3.1) 

where U{pi,p2, q) is the Fourier transform of the interaction Hamiltonian and w are the 
two-component spinor plane-wave solutions to the Dirac equation which is generally 
spin-dependent. The coordinate space representation, r, is the Fourier transform of q 
in the center of mass frame, p = pi = —ph,- After including the relativistic corrections 
to energy and the higher order terms from the Born amplitude, there is a common 
correction to all constituent configurations given by 

AJ/_ = -(i5 + i3)|--^('l:<^#:P + P^y (3.2) 

The primary challenge of computing the leading order correction to the fine structure 
is that it is necessary to use both first and second order perturbation theory. The need 
for second order perturbation theory arises from the fact that the matrix element for 
gaugino exchange between two bound states is parametrically given by 



(1) gaugino _ /ji-ir \>\ _ 



augino I 



U d\— ^\ - a'ml (3.3) 



where the non-relativistic limit of gaugino interaction is described in Sec. 2.2. Due 
to parity, there is no leading order contribution from this interaction^. However, at 
second order 



^^(2) gaugino ^ | | ^gaugino I^Q P ^ «^"^e ^3 



These second order effects cancel off terms from first order perturbation theory that 
would have led to supersymmetry breaking hyperfine interactions. Furthermore, the 
full spectrum of the hydrogen atom contributes (including the continuum states) to 
the second order correction. Schwinger's exact 0(4) symmetric Green's functions were 
used in [1] to get an exact, supersymmetric answer. 



^Dispersive contributions from ~ (f' /nip are further suppressed. 
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Since the dynamics are supersymmetric and the final result must be supersymmet- 
ric, there a is way to side-step the second order calculation and get the exact answer 
from first order perturbation theory. In each of the multiplets of a given energy level, 
there is a state that receives no second order contribution to the energy eigenvalues. 
Therefore, the first order calculation is the final 0{a^) answer. 



Figure 1: The relevant interactions that contribute to the interaction Hamiltonian for the 
first-order exact states: v, <;_, a;_ 

Some of the states that do not receive corrections to their energy eigenvalues from 
second order in perturbation theory are the following states 

\nlm;lv) \nlm]luJ) (3.5) 

where the product between \nlm) ® \S) is performed to go to the total angular mo- 
mentum basis. Only the j = I component of the \lm)\v) state doesn't receive a fine 
structure correction from second order perturbation theory. There are no states at 
second order that can mix these states to because gaugino exchange has the 

following selection rules 

Al = ±1 As = ±1 Aj = (3.6) 

and the states must be of the same parity. These states are representatives of the 
Clifford vacua 

\nlm;lv), \nlm\lq_) G \nl)Vti_i © \nlm;luj^) G \nl)[2 x Qi]. (3.7) 

Therefore, these states represent elements of every independent superfield and by com- 
puting their energy eigenvalues, the entire spectrum is computed. The perturbing 
Hamiltonian has non-vanishing matrix elements between these three states and needs 
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to be diagonalized. Examining the superspin wave functions of these states, only four 
interactions contribute to U{pi,'p2,q) and are illustrated in Fig. 1. After evaluating 
the interactions, there is a 3 x 3 Hamiltonian to diagonalize; however, only |nZm; Iv) 
and |n/m;Z^_) mix. The interaction matrix is as follows 



4 

a jjL 



j^FS Univ. _|_ j^FS Split 



(3.8) 



where 



\nlm] Iv) 
\nlm] /^_) 
\nlm] lujJ) 



FS Univ. 



1 

n 



fx 



(3.9) 



and 



FS Split 



1 



nl 



21 + 1 



2(2/ + 1) 



1(1+1) yz{m) 
^ 







(3.10) 



0/ 



After diagonalizing the interactions and identifying the states with their respective su- 
permultiplets from Eq. 3.7, it is possible to exchange the eigenvalues of orbital angular 
momentum, /, in the expressions for the total angular momentum of the Clifford vac- 
uum, j {e.g. I = j for |n/m; luj_)). The mixing angles between and v are precisely 
the Clebsch- Gordon coefficients for going between the two bases described in Eq. 2.21. 
After performing this substitution, the final fine structure of the supersymmetric hy- 
drogen atom has a remarkably simple (and familiar) expression given by 



nj 



E^y^ + AEf Univ. ^ ^^^^ 



FS Split 



(3.11) 



where the leading order bound state energy is E^^'^ = — a^yu/2n^. The fine structure is 
divided into a j-independent universal term and a term that leads to a fine structure 
splitting 



AE, 



FS Univ. 



4 



[rrip + me] 



A 771FS Split 
nj 



4 



2n^ j + 



(3.12) 



2. As expected, the fine structure 
2.30 when mg/mp — > 0. The next 



The spectrum is schematically illustrated in Fig 
reduces to the non-supersymmetric limit of Eq. 
feature is that the spectrum is independent of the orbital angular momentum, /, as is 
the case in the hydrogen atom. The last pertinent feature is that the j-dependent, fine 
structure splitting receives no supersymmetric contributions. 
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Qs O3 2 X ^ 

Oi 2xfii ^ ^ 

2 



O3 

2 X 1^1 

i 



2 X Oo ^ 

► 

/ = / = 1 1 = 2 

Figure 2: The spectrum of supersymmetric hydrogen up to principal quantum number, 
n = 3. The horizontal horizontal axis denotes the orbital angular momentum, /, beneath 
while above the superspin state, is denoted above. The states are labeled by which 
Clifford they belong to. The dotted lines show the degeneracies that persist after the fine 
structure is incorporated. 



3.1 Supersymmetric Lamb Shift 

The degeneracies of the supersymmetric hydrogen atom at principle quantum number 
n are 2n — 1 two-fold degeneracies. The two-fold degeneracies for the states that are 
built from integer spin Clifford vacua arise from the 0(2)^^ symmetry and unless parity 
or the U{l)ji is broken by some other interactions, these degeneracies will persist to all 
orders in the relativistic expansion as well as to all orders in perturbation theory. 

The two-fold degeneracies for the half-integer Clifford vacua is an accident and it is 
natural to ask at when is the accidental degeneracy lifted. This degeneracy is completely 
analogous to the corresponding degeneracy for the Dirac hydrogen atom and it is well 
known that the Lamb shift arises radiatively. In this section, the correction from 
the Lamb shift is estimated. This article only addresses the logarithmically enhanced 
correction in the infinite proton mass limit. The infinite proton mass limit is particularly 
illuminating because the accidental degeneracy occurs for supermultiplets that have 
valence fermionic electrons. The logarithmically enhanced contribution traces back 
to the collinear singularity of the photon and electron and doesn't exist for scalars 
external legs or photino interactions. Additionally there are fewer interactions because 



E 



^2 y 



2 X i\ 2_ 



2 X On 
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all supersymmetric interactions drop out as discussed in Sec. 2.3. The finite correction 
and corrections away from the infinite proton mass are left for further study. 

In QED there are two one-loop processes contributing to the Lamb shift and the 
IR behavior of the integral is challenging due to the coUinear and soft divergences in 
QED. These divergences are regulated by the discreteness of the hydrogen spectrum. 
Air is the scale where the IR divergences are cut-off. Parametrically Air ~ a^mg and 
detailed calculations give Air ^ lOa^me. 




Figure 3: The SQED contributions to the Lamb shift. Only the upper left diagram has 
a logarithmically enhanced contribution to the splitting. This contribution dominates and 
therefore the dominant contribution to the Lamb shift is identical between QED and SQED. 



All corrections from gaugino exchange and contact interactions from D-terms van- 
ish in the non-relativistic limit as rrip — > oo. The renormalization of supersymmetric 
exchange diagrams can not be parametrically enhanced. The self energy contribu- 
tion from either the photon or photino do not have IR-divergences because the virtual 
particles in the loop are massive. The gaugino corrections to the vertex do not have 
IR-divergences. Therefore, the supersymmetric contributions to the logarithmically 
enhanced portion of the Lamb shift reduce to the QED contributions. Therefore the 
dominant contribution to the Lamb shift is given by [9] 



~ — ^ log ^6io + ■■■ . (3.13) 

The terms that are not logarithmically enhanced will not have the same universal form 
and need to be computed. 

4. Non-Abelian Bound States 

Hydrogen is a useful model for many other bound states and provides intuition about 
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more complicated systems. For instance, the constituent picture of mesons and baryons 
is heavily based on Coulombic bound states of quarks as the fine structure constant 
goes towards unity. This section lays out how the formalism above is used for other 
bound state systems. 

The most straightforward extension is to promote the binding U{1) to an SU{Nc) 
gauge theory with the electron and proton becoming fundamentals and anti-fundamentals 
of the gauge group. When the confinement scale A <C me,mp, the electron-proton 
bound states should be approximately Coulombic. The effective coupling between the 
proton and electron is now given by 

«(^) = C,'-^. (4.1) 

The primary difference between the SQED bound states and SQCD bound states is 
that the gauge coupling runs. The running of the gauge coupling breaks the accidental 
0(4) symmetry of the Coulomb problem. Without the 0(4) symmetry, the (9(a^/i) 
energy eigenvalues are depend upon the orbital angular momentum. The bound state 
spectrum of this theory still has a multi-fold degeneracy arising from different superspin 
configurations ^i, and i^;_|_i) for a given orbital angular momentum, /. This 

remaining degeneracy is lifted by the fine structure. The non-Coulombic nature of the 
running coupling constant will make evaluating fine structure through second order 
perturbation theory difficult. However, the method of identifying the states that do 
not have corrections to their energy eigenvalue from second order perturbation theory 
is still applicable. 

A related example to the previous case is to keep the proton and electron heavy 
and add in Nf light spectator quarks that slow the running of a(/i). In the limit 
Nf ~ 3Nc the theory enters a weakly coupled fixed point in the IR and the leading 
order potential becomes approximately Coulombic and the 0(4) symmetry re-emerges. 
The arguments from Sec. 2.2 indicate that the fine structure degeneracy will emerge 
in the infinite proton mass limit. The Lamb shift is more subtle in this case because 
the gluons of the SU{Nc) gauge theory are colored. The logarithmically enhanced 
contribution is related to emitting a gluon; however emitting a gluon from the binding, 
color-singlet state will transition the bound state into a color-adjoint configuration that 
doesn't bind. The virtual gluon has energy of the order of a Rydberg and the energy- 
time uncertainty softens this effect. Nevertheless, this discussion indicates that the 
Lamb degeneracy is potentially more interesting for these non-Abelian gauge theories. 

4.1 Heavy-Light Bound States 

The final example of a non-Abelian bound state that this article considers is a more 
subtle situation when there is one light quark (the electron) and one heavy quark (the 
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proton), relative to the scale of confinement. Non-Abelian super- Yang Mills theories 
with Np < Nc do not have a vacuum in the absence of supersymmetry breaking. 
Instead, the hght quarks develop an Affleck-Dine-Seiberg superpotential that is singular 
at the origin 

/ A^ \ ^ 

^ADs - (^J (4.2) 

where E and E'^ are the light-flavored quarks. In the presence of soft supersymmetry 
breaking, soft masses for the squarks will stabilize the run-away potential. The location 
of the vacuum is roughly 

/ A \ 2{iVc + l) 

(g) = (g^) ~ A ~ my. (4.3) 

V"^soft/ 

These scalars can only break SU{Nc) SU{Nc — 1) and it gives mass to vector 
superfields transforming as □ + □ -|- 1 under the SU{Nc — 1). The electrons get eaten 
by the super-Higgs mechanism and absorbed into the massive vectors. This stabilizes 
the gives mass to the vectors, my, above the confinement scale so long as A > rrisoft- In 
this is the case, the massive vectors are quasi- weakly coupled and the resulting spectrum 
of the theory is qualitatively different than the those of supersymmetric hydrogen. 

The heavy proton, after the scalar electron acquires a vacuum expectation value, 
decomposes into 

P ^ 1 + □, 1 + □. (4.4) 

Since the SU{Nc — 1) has no light flavors, it confines at low energies and the colored 
proton must be screened. The proton can only be screened by an anti-proton or one 
of the massive vectors (that have eaten the electrons). Because the vectors are heavier 
than the confinement scale, the proton-vector bound states are quasi-Coulombically 
bound. The procedure from Sec. 2.1 can be performed again, but with \S) now given 
by the product of a hypermultiplet and a massive vector field: 

= (P © P'^)p-p>„„+p-,,/2 ® V^p=p„„-p-.,/2 (4.5) 
This product can be decomposed into 

= rii massive, complex spin | supermultiplet 

= 2 X f2i massive, complex vector supermultiplet (4-6) 
= ilo massive, complex hypermultiplet 
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The 0(2)r symmetry is still present in the theory and causes the doubling of the vector 
supermultiplet. Since the ground state will have no orbital angular momentum, \S) are 
equivalent to the ground state degeneracy of the system. The fine structure will split 
these multiplets. 

In addition to the supersymmetric corrections to the spectrum, supersymmetry 
breaking effects arise from the soft masses of the theory; however, these can be made 
arbitrarily small in the limit msoft/A — > holding mv/rrip fixed. 

5. Discussion 

This article calculated the spectrum of supersymmetric hydrogen and yielded a spec- 
trum that is nearly identical to that of the Klein- Gordon and Dirac equation even for 
finite proton masses. The surprisingly mild effect of the plethora of supersymmetric in- 
teractions is closely related to the non-relativistic limit. The fine structure degeneracy 
is still present in the supersymmetric hydrogen atom and this degeneracy is lifted by the 
Lamb shift. The logarithmically enhanced contribution to the Lamb shift is the same as 
in QED, but the finite contributions are sensitive to supersymmetric interactions. The 
Johnson-Lippmann operator guarantees the / independence of the energy eigenvalues in 
the infinite proton mass limit for both supersymmetric and non-supersymmetric QED. 
For a finite proton mass in non-supersymmetric QED there is no degeneracies of the fine 
structure and therefore no Johnshon-Lippmann operator. However, in supersymmetric 
QED, the residual degeneracy of the fine structure indicates that a generalization of 
the Johnson-Lippmann operator should exist even with a finite proton mass and this 
operator has not been found. The methods presented in this article are widely appli- 
cable to other supersymmetric bound states and additional applications were briefiy 
discussed in Sec. 4. 

In addition to studying the spectroscopy of atoms, supersymmetry also constrains 
the interactions of the bound states. The interactions of the ground state with an 
external fields (such as the vector superfield) should be described in terms of an effective 
field theory. Studying the interactions may giver further insight into the dynamics of 
the bound states. Ultimately, some variant of these bound states might be dark matter. 
The astrophysical anomalies currently observed hint that dark matter may have more 
structure than typically assumed, see [10] for a recent description. The tools presented 
in this article may help construct models to further explore these possibilities. The 
theory studied in Sec. 4.1 is the supersymmetric generalization of composite inelastic 
dark matter [4,6]. The discussion in Sec. 4.1 indicates that there is a much richer 
structure than the simple generalization of mesons to supersymmetric mesons because 
of the vacuum structure of the theory. 
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If this supersymmetric gauge theory is coupled to the Standard Model, supersym- 
metry breaking will be communicated to the spectrum at some level. At this stage, it is 
not possible to cherry pick convenient states to compute the effects of supersymmetry 
breaking; however, assuming that supersymmetry breaking is a modest effect, it is only 
necessary to use first order in the supersymmetry breaking interactions. The leading or- 
der splitting is induced by the scalar soft masses. This perturbation is straight-forward. 
The more challenging measure is the two-fold degeneracy of states in the \S) = 
configuration will only be lifted through parity breaking or i?-breaking. These effects 
are model dependent and a specific example is considered in [7]. 
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